Diffusion and viscosity in a supercooled polydisperse system 
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We have carried out extensive molecular dynamics simulations of a supercooled polydisperse 
Lennard- Jones liquid with large variations in temperature at a fixed pressure. The particles in the 
system are considered to be polydisperse both in size and mass. The temperature dependence of 
the dynamical properties such as the viscosity (77) and the self-diffusion coefficients (Di) of different 
size particles is studied. Both viscosity and diffusion coefficients show super-Arrhenius temperature 
dependence and fit well to the well-known Vogel-Fulcher-Tammann (VFT) equation. Within the 
temperature range investigated, the value of the Angell's fragility parameter (D ~ f.4) classifies the 
present system into a strongly fragile liquid. The critical temperature for diffusion (T®*) increases 
with the size of the particles. The critical temperature for viscosity (T%) is larger than that for 
the diffusion and a sizeable deviations appear for the smaller size particles implying a decoupling 
of translational diffusion from viscosity in deeply supercooled liquid. Indeed, the diffusion shows 
markedly non-Stokesian behavior at low temperatures where a highly nonlinear dependence on size 
is observed. An inspection of the trajectories of the particles shows that at low temperatures the 
motions of both the smallest and largest size particles are discontinuous (jump-type). However, 
the crossover from continuous Brownian to large length hopping motion takes place at shorter time 
scales for the smaller size particles. 



I. INTRODUCTION 



The rapid cooling of a liquid below its freezing temper- 
ature transforms it into a long-lived metastable amor- 
phous solid or glass Understanding the dynamics 
of the system near the glass transition is an intense 
field of research since the last few decades. There have 
been many experimental 3|4 ,|5|||(7| as well as simula- 
tion studies [||Jl^,|ll|,|l|,y,y| which focus on the dy- 
namics of dense supercooled liquids well above the glass 
transition and also near the glass transition temperature 



(T g ). The basic aim of all these studies was to charac- 
terise quantitatively the observed very complex dynamics 
of the system as it approaches the glass transition from 
the above. Close to the glass transition, the shear vis- 
cosity (77) and the microscopic structural relaxation time 
(t) of the so-called fragile glass-forming liquids show di- 
vergence with a strongly non-Arrhenius temperature de- 
pendence Q. This divergence is often well represented 
by the Vogel-Fulcher-Tammann (VFT) equation 



fj(T) = A v exp[E n /(T-T2)} 



(1) 



where and E v are temperature- independent constants 
and T 7 J (< T g ) is the temperature at which r\ diverges. 
Note that at low temperatures the increasingly slow dy- 
namics of the so-called fragile liquids is simultaneously 
manifested in the stretched exponential decay of the 
stress correlation function (with a strongly temperature 



dependent stretching parameter) |]15| , ^6| . The VFT de- 
pendence (Eq. 1) is thus accompanied by the strong non- 
exponential relaxation observed near the glass transition 
temperature. 

The dramatic slow down of the dynamics near the 
glass transition is not well understood and still remains 
the most challenging problem in the physics of glasses. 
Several theories have been proposed to understand the 
anomalous relaxation dynamics of deeply supercooled liq- 
uids. Although the ideal version of the non-linear mode 
coupling theory (MCT) |l7| gives a microscopic picture 
of this slowing down, it predicts a structural arrest, i.e., 
a transition from ergodic to nonergodic behavior, at a 
critical temperature T c , well above the laboratory glass 
transition temperature T g . Near T c , the importance of 
the influence of potential energy landscape on the relax- 
ation processes has now widely been accepted Jl8|,[l9| and 
the strongly correlated jump motion is observed to be the 
dominant mode for mass transport p3| , p0| , pl 22 1, which 
is not included in ideal MCT. 

Another important characteristic feature in the dy- 
namics of deeply supercooled liquids is the decoupling 
between translation diffusion and the shear viscosity of 
the medium |2^,^J,^5|. At high temperature, over a wide 
range of liquid states, the translational diffusion is in- 
versely proportional to viscosity, in accordance to the 
Stokes-Einstein (SE) relation given by 



D T = 



k B T 
CirrjR 



(2) 
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where R is the spherical radius of the diffusing par- 
ticle and C is a numerical constant that depends on 
the hydrodynamic boundary condition. However, sev- 
eral recent e xperimental |26|j27|j2^] and simulation studies 
i^,|33[| on strongly supercooled 'fragile' glass 
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forming liquids have shown significant deviations from 
the SE relation. As the temperature is lowered toward 
T g , it is found that the translational diffusion is larger 
than the value predicted by the SE relation and in some 
cases even two to three orders of magnitude larger Q. 
The enhanced diffusion at low temperatures is sometimes 
explained in terms of a power law behavior Dt oc r\~ a 
with a < 1 [28 3C(]. Both the experiments and simulation 



studies have recently been evidenced the enhancement of 
the translational diffusion coefficient is due to the spa- 
tially hete rogen eous dynamics in deeply supercooled liq- 
uids |,0|H|,||. 

Computer simulations have played a key role in aug- 
menting our understanding of various aspects of the dy- 
namics of the supercooled liquids from a microscopic 
viewpoint. Unfortunately, the simple one-component 
systems such as soft or hard spheres or Lennard-Jones 
systems crystallize rapidly on lowering the temperature 
below the melting point (T m ) and, therefore, cannot be 
utilized as a model for studying the complex dynami- 
cal behavior near the glass transition temperature. A 
natural way to avoid crystallization is to use binary mix- 
tures of atoms with different diameters. A large number 
of molecular dynamics (MD) simulations have recently 
been carried out in supercooled model binary mixtures 
near the glass transition as well as below the glass tran- 
sition temperature g|@||||,0,P,||||,||@ . 

However, one is often interested in the consequences 
of the disorder introduced by the dissimilarity of the 
particles. Synthetic colloids, by their very nature, fre- 
quently exhibit considerable size polydispersity fl37| , |3"8| | . 
Polydispersity is also common in industrially produced 
polymers which always contain macromolecules with a 
rang of chain length. Colloidal particles are an excellent 
model of hard spheres and perhaps the simplest possible 
experimental system of interacting particles to study the 
glas s transition. Several experiments [ p7[ and simulations 
[p9| , ^0| have shown that the crystal phase of the colloidal 
systems can exist as a thcrmodynamically stable phase 
only for polydispersities (standard deviation of the size 
distribution divided by the mean) less than a 'terminal' 
value, in the range of 0.05-0.15. 

Interestingly, recent experiments on colloidal super- 
cooled fluids and colloidal glasses allowed one to obtain 
the information on the microscopic details of the dynam- 
ics of the individual particles || . These experiments have 
shown the presence of dynamic heterogeneity in deeply 
supercooled colloidal systems. The motion of the rela- 
tively fast moving particles is found to be highly corre- 
lated and form connected clusters, whose size increases 
as one approach the glass transition. More recently, Sear 
pf[ has carried out MD simulation of a dense polydis- 
perse hard sphere fluid to study the effect of polydis- 



persity on the slow dynamics. The simulation results 
also show the clustering of the fast-moving particles in 
agreement with the experiments, though the dynamics 
appears to be less heterogeneous. The heterogeneous na- 
ture of the dynamics has also been observed in the Monte 
Carlo simulation study of polydisperse hard spheres close 
to the glass transition (42] . 

It is worth mentioning that besides the short-range 
hard-core interaction, addition of a nonadsorbing solu- 
ble polymer in the stable colloidal suspension gives rise 
to a weak, long-range attraction between the colloidal 
particles by means of the depletion interaction |}8| . Re- 
cently, the consequences of this attractive interaction on 
the glass transition are nicely explained in a combined 
experimental, theoretical, and simulation study by Pham 
et.al. |fj| Interestingly, with increase in the strength of 
the short-range attractive interaction, two qualitatively 
different glassy states are found with a reentrant glass 
transition line. 

The size distribution in real colloids generally leads 
to a distribution in mass of the particles. The impor- 
tance of the mass polydispersity on the dynamics of a 
realistic system having size polydispersity has recently 
been analyzed in molecular dynamics simulation study of 
a Lennard-Jones (LJ) polydisperse fluid near the triple 
point of the corresponding monodisperse L J system p4[ . 
Polydispersity is commonly found in many systems of in- 
dustrial applications and to mimic the interparticle inter- 
actions, the Lennard-Jones potential generally serve as a 
good starting model. Thus, it will be a general interest 
to study the impact of polydispersity on the dynamics of 
a deeply supercooled polydisperse fluid, where particles 
interacting via the LJ potential. More importantly, this 
will enable us to compare the properties of the system 
with model binary LJ mixtures whose dynamics near the 
glass transition has been studied extensively in simula- 
tions. 

In this study, we have performed extensive molecu- 
lar dynamics simulations of a system of polydisperse LJ 
spheres with continuous range of diameters and mass. 
The temperature dependence of the dynamic properties 
such as viscosity (rj) and the self-diffusion coefficients 
(Di) for the different size particles is studied by vary- 
ing the temperature (T) over a large range at a con- 
stant high pressure (P). Both the viscosity and diffusion 
show super-Arrhenius temperature dependence and the 
calculated value of the fragility parameter (D) show that 
the present system is more fragile than the well-known 
Kob- Andersen binary mixture. The critical temperature 
obtained from the VFT fit to the diffusions (T^ 4 ) show 
strong dependency on the radius (i?) of the particles. In 
addition, the critical temperature obtained from the VFT 
fit to the viscosity (T%) is much higher when compared 
to those of diffusion coefficients, where the deviation is 
largest for the smaller size particles. This clearly reflects 
the deviation from the Stokesian diffusion in the prox- 
imity of the glass transition temperature. Most interest- 
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ingly, at lower temperatures the diffusion shows a highly 
nonlinear size dependence when plotted against the in- 
verse of the radius (R) of the particles. The reason for 
the breakdown of the Stokes-Einstein (SE) relation can 
be analyzed from the trajectories of the particles. We 
find that at low temperature the hopping processes be- 
ing the primary mode of particle diffusion for both the 
smaller and bigger size particles. 

The organization of the rest of the paper is as follows. 
In the next section, we describe in detail the system stud- 
ied here and the details of the simulations. The simu- 
lation results are analyzed and discussed in section III. 
Finally, we end with some concluding remarks in section 
IV. 



II. SYSTEM AND SIMULATION DETAILS 

We have performed a series of equilibrium isothcrmal- 
isobaric ensemble (N P T) molecular dynamics simula- 
tions in three dimensions of a system of N = 256 par- 
ticles of mean radius a and mass fh with polydispersity 
in both size and mass. The interaction between any two 
particles is modeled by means of shifted force Lennard- 
Jones (LJ) pair potential, where the standard LJ is given 
by §1 
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where i and j denotes two different particles and a^j = 
((Tj +r7j)/2 where <7j, (jj are the diameters of the particles 
i and j, respectively. In the shifted-force potential both 
the potential and force are continuous at a cutoff radius 
r c and we choose a value of r c = 2.5a\ 

The polydispersity in size is introduced by random 
sampling from the Gaussian distribution of particle di- 
ameters a 1441 



P(a) = 



<5V2tt 



exp 




(4) 



where 8 is the width of the distribution. The standard 
deviation (5) of the distribution divided by its mean a 
gives a dimensionless parameter, the polydispersity in- 
dex s = 6/ a. The simulations are carried out here for 
a fixed value of the polydispersity index s — 0.1. The 
masses of the particles are varied accordingly and we as- 
sume that the mass of a particle i is scaled by its diameter 
as rtii — fh(ai/a) 3 . The LJ energy parameter is as- 
sumed to be same for all particle pairs and denoted as e. 
All the quantities in this study are given in reduced units, 
that is, length in units of a, temperature T in units of 
e/fe_B, pressure P in units of e/a 3 , and the unit of time is 
t = y/ ma 2 j t. Note that if one assumes the argon units 
then r = 2.2 ps. 



All simulations in the NPT ensemble were performed 
using the Nose-Hoover- Andersen method jl6| , where the 
external reduced temperature (T*) is varied over a large 
range from 1.3 to 0.67 keeping the external reduced pres- 
sure (P*) fixed at 10.0. Throughout the course of the 
simulations, the barostat and system's degrees of freedom 
are coupled to an independent Nose-Hoover chain Ji?]] 
(NHC) of thermostats, each of length 5. The extended 
system equations of motion are integrated using the re- 
versible integrator method |)48|| . The higher order mul- 
tiple time step method has been employed in the NHC 
evolution operator which lead to stable energy conserva- 
tion for non-Hamiltonian dynamical systems p9fl . The 
extended system time scale parameter used in the calcu- 
lations was taken to be 0.93 for T* > 1.0 and 1.16 for 
T* < 1.0 for both the barostat and thermostats. 

A time step of 0.001 is employed for T* > 1.0 and 0.002 
for T* < 1.0. The equilibration and data collection steps 
are also varied accordingly depending upon the tempera- 
ture of the system. For T* > 1.0, the equilibration steps 
are varied from 5 x 10 5 to 10 6 and the data collection steps 
are 10 6 , whereas for T* < 1.0, the equilibration steps are 
varied from 5.0 x 10 5 to 2.0 x 10 6 and the data collec- 
tion steps from 10 6 to 2.5 x 10 7 . At each temperature, 
all the dynamic quantities are averaged over five indepen- 
dent runs. Diffusion coefficients (Di) for the different size 
particles are calculated from the slope of the correspond- 
ing mean square displacements (MSD) in the diffusive 
limit and viscosity is calculated from the auto-correlation 
of the off-diagonal components of the microscopic stress 
tensor, via the standard Green-Kubo formula |30]]. As 
the system is isotropic, we have taken an average over 
three different off-diagonal stress correlations for each of 
the five data set. 



III. RESULTS AND DISCUSSION 

In order to make sure that there is no crystallization, 
we have calculated the radial distribution functions g(r) 
which describes the average structure of the fluid. The 
radial distribution function calculated for T* — 0.67 is 
shown in figure 1. The decay of correlations with in- 
crease in distance clearly reflects the absence of any long 
range order, a characteristic feature of the fluid. 

The plot of ln(77) as a function of inverse of temperature 
(1/T*) in figure 2(a) clearly shows a super- Arrhenius be- 
havior of the viscosity. In figure 2(b) we show a VFT 
fit to the viscosity (Eq. 1) by plotting ln(^) against 
(1/(T-T^)) where TJj is equal to 0.57. As in other fragile 
liquids, it shows that the divergence of viscosity is quite 
well described by the VFT equation. From the fitting we 
obtain the values of A v and E v as 2.0 and 0.81, respec- 
tively. We have also calculated the fragility parameter (D 
= E v /T£) as defined by Angell Jsljl . Using the values of 
the fitting parameters (E^ and T%) obtained within the 
temperature range investigated, its value is ~ 1.42. This 
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classifies the present system into a strongly fragile liquid. 
Thus, the dense random packing of unequal size particles 
makes the present system more fragile when compared 
with a recent simulation study on Kob- Andersen binary 
mixture (D w 2.45) 

In a polydisperse system, all the particles arc unequal 
in size, so their diffusion coefficients will also differ. We 
categorize the particles into different subsets where parti- 
cles of diameters within 0.055- are assumed to be members 
of the same subset. For the polydispersity index s = 0.1, 
we find that the minimum and maximum diameter of the 
particles are 0.75a and 1.25a, respectively. Thus subsets 
of particles with diameters in the ranges 0.75 to 0.8a and 
1.2 to 1.25(7 corresponds to smallest and largest spheres, 
respectively. The diffusion coefficients for different sub- 
sets of particles are calculated at each temperature. It is 
well-known that in deeply supercooled liquids the diffu- 
sion coefficient shows non-Arrhenius temperature depen- 
dence and can be fitted by a VFT law 

Di(T) = A Di exp[-E D J{T-T°% (5) 

where the index i stands for the different subsets of par- 
ticles. T^ >i is the critical temperatures for i-th species at 
which the diffusion coefficients (Di) vanishes. The dif- 
fusion coefficients for each subsets of particles have been 
fitted to the above equation and we show the plot of 
ln(A) against 1/(T - T® % ) , for the smallest (D x ) and 
largest spheres (£>io) in figures 3(a) and 3(b), respec- 
tively, where 1 and T Dl0 are 0.46 and 0.5. Thus the 
critical temperatures (T^*) depends on the size of the 
particles and increases with size of the particles Ri (Ri 
is the mean radius of the i-th subset). This is shown in 
figure 4. For the largest spheres, the critical temperature 
(T® 10 ) is smaller than the corresponding critical temper- 
ature obtained from VFT fit to the viscosity (Tg = 0.57). 
This clearly signifies that near the glass transition the 
diffusion is partly decoupled from the viscosity and for 
smaller particles the degree of decoupling is more. The 
smaller particles remain mobile even when bigger parti- 
cles are almost frozen. 

In figure 5 we plot the diffusion constants (Di) against 
1/Ri at the lowest temperature of T* = 0.67 and 
compared with the well-known hydrodynamic Stokes- 
Einstein (SE) relation (Eq. 2 with C = 6, the stick 
boundary condition) . It clearly shows the markedly non- 
Stokesian behavior of the diffusion at low temperatures. 
Interestingly, the fitting to the simulated data points 
show a highly nonlinear size dependence of the diffusion. 
This is a clear evidence that the breakdown of SE law 
is more severe for the smaller size particles. In order to 
get an estimate of the degree of decoupling (between dif- 
fusion and viscosity) for the smallest size particles, we 
have fit the inverse diffusion coefficient (1/D\) versus 
rj/T. While at high T, it asymptotically satisfies the SE 
relation (slope is 1), the fit to the low temperature data 
gives the slope a « 0.5 (that is, diffusion shows the power 
law behavior D\ oc 7y -0 - 5 ). It should be noted that the 



dynamics of a polydisperse liquid is more heterogeneous 
than a monodisperse or bidisperse system due to the dif- 
ferent time scales involved for different sizes and masses 
of the particles Q^^] . The smaller particles are on aver- 
age faster than others over all time scales. This becomes 
more prominent at lower temperatures where the relax- 
ation time of the system is very high. At low tempera- 
ture, the observed nonlinear dependence of diffusion on 
size is related to the increase in dynamic heterogeneity 
in a polydisperse system. 

A more detailed analysis of the diffusion can be ob- 
tained from a closer examination of the self-part of the 
van Hove correlation function G s (r, t) . This gives the 
distribution of the displacements (r) of a particle in a 
time interval t. We calculate G s (r, t) for the smallest (a 
= 0.75 to 0.8a) and largest (a = 1.2 to 1.25a) particles 
for different time intervals at T* = 0.67, the lowest tem- 
perature investigated. Figures 6(a) and 6(b) shows the 
correlations for the smallest and largest spheres, respec- 
tively. For the smallest particles, a gradual development 
of a well-defined second peak at r ~ 1.0a is clearly visi- 
ble with increase in time (figure 6(a)). However, for the 
largest particles, the distribution becomes bimodal at rel- 
atively longer time scales (figure 6(b)). The occurrence of 
the secondary peak, observed also in other model binary 
mixtures at low temperatures 0,^3^51, is an evidence 
of the jump motion in the dynamics of the particles. 

To characterize the single particle dynamics further, 
we have evaluated the self-intermediate scattering func- 
tion F s (k,t), the spatial Fourier transform of G s (r,t), for 
different subsets of the particles for a fixed value of the 
reduced wave number ka ~ 2tt at T* = 0.67. The long 
time decay of F s (k,t) is well fitted by the Kohlrausch- 
Williams- Watts (KWW) stretched exponential form 

Fi(k,t) = exp^-^ (6) 

where Tj and [3i are the relaxation time and the stretching 
exponent of the i-th subset. We find that both r, and Pi 
increases with an increase in the size of the particles, as 
has been observed earlier by other authors in binary mix- 
tures [0. The F s (k, t) calculated for the smallest (subset 
1) and largest (subset 10) particles along with the KWW 
fits is shown in figures 7(a) and 7(b), respectively. Note 
that we fit the functions [F l s (k, t - t )/F l s (k, t )] (t > t a ), 
to the KWW form to quantify their long time behavior 
||. For the smallest particles, the values of the fitting 
parameters are found to be t\ ~ 242 and P\ ~ 0.49, 
whereas for the largest particles they are tio — 717 and 
/3io — 0.64. The enhanced stretching (Pi < pw) at long 
times is due to the more heterogeneity probed by the 
smaller size particles than that by the larger size parti- 
cles during the time scale of decay of their F s (k, t) |5g ]. 
In order to determine the extent of the jump-type mo- 
tion more clearly we follow the trajectory of the individ- 
ual particles. A close inspection of the simulated trajec- 
tory of the smallest and largest particles reveals several 
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interesting features. Figure 8(a) display the projections 
onto an x-y plane of the trajectory of a typical smallest 
size particle over a time interval At = 500r and figure 
8(b) shows the paths followed by a largest size parti- 
cle over time interval At = 2000r, both at T* = 0.67. 
At this temperature the dynamics is dominated by 'hop- 
ping'; particles remain trapped in transient cages created 
by the surrounding particles for quite some time and then 
moves significant distance (approximately one interpar- 
ticle distance) by making a jump to a new cage. For 
the larger size particles, the jump motion begins to take 
place at relatively longer time scales (compare also fig- 
ures 6(a) and 6(b)). Thus, it clearly elucidates that in 
a system with particles of all different sizes and masses, 
the hopping is the dominant diffusive mode both for the 
smaller and bigger size particles. It is the frequent hop- 
ping in case of smaller size particles that leads to the 
severe breakdown of SE relation. 



IV. CONCLUSIONS 

In summary, we have presented the results of large 
scale computer simulations for a supercooled polydis- 
perse system with large variations in temperature at a 
fixed high pressure. The interparticle interaction is rep- 
resented by the standard Lennard- Jones (LJ) potential 
and the particles in the system are considered to be 
polydisperse both in size and mass. Characteristic of 
a fragile glass former, the super-Arrhenius temperature 
dependence is observed for the viscosity and also for the 
self-diffusion coefficients of different size particles. Fur- 
thermore, within the temperature range investigated, the 
value obtained for the Angell's fragility parameter (D 
w 1.4) establish the present system as a strongly fragile 
liquid. 

In a dense polydisperse system, a wide range of time 
scales are involved due to the different size and mass of 
the particles. Thus, upon lowering the temperature, the 
dynamics is increasingly more heterogeneous compared 
to a monodisperse or bidisperse system. The increase of 
critical temperature for diffusion (T®*) with the size of 
the particles suggests that the dynamics is indeed hetero- 
geneous. In addition, the critical temperature for viscos- 
ity (T%) is found to be larger than that for the diffusion, 
indicating the decoupling of translational diffusion from 
the viscosity commonly observed in deeply supercooled 
liquids. Interestingly, a marked deviation from the Stoke- 
sian diffusion is observed where the dependence on size 
of the particles is highly nonlinear. 

The analysis of the self-part of the van Hove correla- 
tion functions G s (r, t) showed a clear signature of the sin- 
gle particle 'hopping' (indicated by the second neighbor 
peak) in the dynamics of both the smallest and largest 
size particles at low temperatures. However, for the 
larger particles, the hopping processes are found to occur 
at relatively longer times. The relevance of the hopping 



processes at low temperatures is further investigated in 
detail by following the trajectory of the individual par- 
ticles. Crossover from continuous Brownian to hopping 
motion takes place at shorter time scales for the smaller 
size particles than that for the larger size particles. 

In the present system the size of all the particles are 
different. It would be interesting to sec whether the jump 
motion executed by the individual particles occurs over 
a single energy barrier or it takes place via a number of 
"intermediate" inherent structures in the potential en- 
ergy landscape. A recent molecular dynamics simula- 
tions on LJ binary mixture |lj] have shown that such 
a transition does not correspond to transitions of the 
system over single energy barriers. In addition, in a 
deeply supercooled liquid the jump motions are associ- 
ated with strong nearest-neighbour correlations, in which 
several neighboring atoms jump at successive close times 
po| , pTlj35| . It is to be noted that similar correlations have 
been observed here also. Recently, a computer simula- 
tion study of a deeply supercooled binary mixture ]2^] 
has shown that the local anisotropy in the stress is re- 
sponsible (at least partly) for the particle hopping. How- 
ever, the molecular origin of the jump motions observed 
here (highly disordered system) is not clear and we are 
presently pursuing this problem. 
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Figure Captions 

Figure 1. The radial distribution function g(r) of 
the system at T* — 0.67, the lowest temperature investi- 
gated. 

Figure 2. Temperature dependence of the shear vis- 
cosity (-q). (a) ln(?7) as a function of inverse of temper- 
ature (1/T*). The simulating values given by the solid 
circles show super- Arrhenius behavior. The dashed line 
gives a guideline to the Arrhenius behavior, (b) 111(77) 
against — TJ). The solid circles are again repre- 

sents the simulation results and the VFT fit is shown by 
the solid line. T% is found to be 0.57. The slope (E v ) 
and intercept (In A v ) obtained from the fit are 0.81 and 
0.69, respectively. For details, see the text. 

Figure 3. Temperature dependence of the diffusion 
coefficients for the smallest (D\) and largest (Dio) par- 
ticles, (a) ln(£>i) is plotted against 1/(T - T® 1 ). (b) 
ln(Dio) against 1/(T -T^°). Solid circles are the sim- 
ulation results and the solid lines are the VFT fit. The 
critical temperatures T® 1 and T^ 10 obtained from the 
VFT fit are 0.46 and 0.50, respectively. 

Figure 4. The critical temperature (T^) obtained 
from the VFT fit to the different subsets of particles as a 
function of the mean radius (Ri) of the subsets (in units 
of er) . Note that T^ >i increases with the size of the par- 
ticles. 

Figure 5. The diffusion coefficients (Di) as a func- 
tion of 1/Ri at T* = 0.67. The dashed line repre- 
sents the Stokes-Einstein relation (Eq. 2) with the stick 
boundary condition C = 6. The viscosity (77) value is 
taken from the present simulations. Solid circles are the 
simulated values and the solid line is the cubic poly- 



nomial fit in 1/Ri. The fit parameters are as follows: 
D, = 0.0011-0.00132(l/i?,)+0.000442(l/i? I ) 2 - It clearly 
shows a highly nonlinear size dependence and a marked 
deviation from the Stokesian behavior of the diffusion. 

Figure 6. The van Hove self-correlation function 
G s (r, t) as a function of the particle displacements r (in 
units of a) at T* — 0.67 for different values of time t (in 
units of t = yj ma 2 / e — 2.2 ps for argon units), (a) For 
the smallest size particles (subset 1). The occurrence of 
the second peak at r w I.Oct indicates the single particle 
hopping, (b) For the largest size particles (subset 10). 
Here also a second peak corresponds to single particle 
hopping develops but at relatively longer times. 

Figure 7. The self-intermediate scattering function 
F s (k,t) for T* — 0.67 is shown with a shift in the time 
origin to t — 1.0, and normalized to the value at t Q , for 
a fixed value of ka = 2tt. This transformation is a conve- 
nient way to eliminate the Gaussian dependence at short 
time |]]. (a) For the smallest size particles (subset 1). 
(b) For the largest size particles (subset 10). Open cir- 
cles represents the simulation results and the solid lines 
are the stretched exponential fit (Eq. 6). The time con- 
stants (n and tio) and the exponents (/?i and /3io) ob- 
tained from the fits are T\ ~ 242, (3i ~ 0.49, r 10 — 717, 
and P w ~ 0.64. 

Figure 8. (a) Projections into x-y plane of the tra- 
jectory of a typical smallest size particle over a time in- 
terval t = 500r. (b) Projections into x-y plane of the 
trajectory of a typical largest size particle over a time 
interval t = 2000r. Note that the time (t) is scaled by 
t = y/ ma 2 1 ' t it is 2.2 ps if argon units are assumed. For 
detailed discussion, see the text. 
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